FINITE-GAP POTENTIAL, HEUN'S DIFFERENTIAL EQUATION 

AND WKB ANALYSIS 



KOUICHI TAKEMURA 

Abstract. We review several results on the finite-gap potential and Heun's dif- 
ferential equation, and we discuss relationships among the finite-gap potential, the 
WKB analysis and Heun's differential equation. 



1. Introduction 

Heun's equation (Heun's differential equation) is a linear differential equation of 
second order given by 

(11) | f 7 | - | e ^ dy | ~ q y-Q 

dz 2 \z z — 1 z — t J dz z(z — l)(z — t) 

with the condition j + 8 + e = a + (3+l [12]. It has four singularities {0, 1, t, oo} and 
they are all regular. Heun's equation is known to be a standard form of the second- 
order Fuchsian differential equation with four singularities. The parameter q is not 
determined by the local monodromy, and is called an accessory parameter. Heun's 
differential equation frequently appears in Physics, i.e. black hole (general relativity, 
Kerr's solution), crystalline materials [13] . fluid dynamics [3], quantum mechanics 
(Inozemtsev model |llj ) and so on. 

A standard form of the second-order Fuchsian differential equation with three sin- 
gularities is given by the hypergeometric differential equation 

(1.2) z(l - z)^+ ( T - (a + (3 + l)z) %-<*Pv = 0. 

It has three regular singularities at {0, 1, oo}. Global properties of solutions and the 
monodromy are known for the hypergeometric differential equation. In particular an 
integral representation of a solution is given by 

(1.3) F(a, ft 7; z) = — ^ r f s«-\l - S y—\l - 8z)-fd8. 

r(a)r(7-a) Jo 

The connection matrix among local solutions at z = and ones at z = 1 is written 
in terms of the gamma function. The condition for existence of non-zero function 
holomorphic at z = and z = 1 is known and the holomorphic function is given by 
the Jacobi polynomial. 

For investigating global properties for Heun's differential equation, we can apply the 
method of finite-gap integration. For this purpose, we recall an elliptic representation 
of Heun's differential equation. Let p(x) be the Weierstrass doubly-periodic function 
with periods (2o>i, 2u 3 ). Set u> = 0, u 2 = — uj\ — u> 3 , = p(u)i) (i = 1, 2, 3), 

(1.4) z= fM^l t ( = fl^£l, = - 

e 2 - ei &2- e 1 
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then Heun's differential equation (Eq. fll.ip ) is transformed to 
(1.5) {-^ + l $i + l)fp(x + Ui) - E^j f(x) = 0, 

where 
(1.6) 

l = (3-a-l/2, / 1 = - 7 + l/2, l 2 = -6 + l/2, h = -e + l/2, 
E=(e 2 - e x ){-Aq + {-{a - (3) 2 + 2 7 2 + 67c + 2e 2 - 4 7 - 4e - 5 2 + 25 + l)/3 
+ (-(a - p) 2 + 2 7 2 + 675 + 25 2 - 4 7 - A5 - e 2 + 2e + l)t/3). 

If Zi = l 2 — h — (7 = 6 = e = 1/2), then the differential equation is called Lame's 
equation. 

It is known that, if Iq, h, l 2 , h € 2, then the function Y^=o h(k + l)p(^ + ^>i) is an 
algebro-geometric finite-gap potential, and is called the Treibich-Verdier potential. In 
section [2} we recall the definitions on the finite-gap potential and review the properties 
and the examples including the Treibich-Verdier potential. The monodromy of Heun's 
differential equation is investigated in connection with the finite-gap property, and 
we will obtain formulae related with the finite-gap property in section [31 In section 
[H we discuss relationships among the finite-gap potential, the WKB analysis and 
Heun's differential equation. In particular, we provide another approach for results 
by Borcea and Shapiro [2] on root asymptotics of spectral polynomials for the Lame 
operator. 

2. Finite-gap potential 

We recall definitions of the finite-gap potential and the algebro-geometric finite-gap 
potential. 

Definition 1. Let q(x) be a periodic, smooth, real function, H be the operator 
—d 2 /dx 2 + q(x), and the set o- b (H) be defined as follows: 

(2.1) E G o~b(H) 45 Every solution to (H — E)f(x) = is bounded on x G M. 
If the closure of the set o~b{H) can be written as 

(2.2) a b {H) = [E , E x ] U [E 2 , E 3 ) U • • • U [E 2g , 00), 

where E < E x < ■ ■ ■ < E 2g , i.e., the number of bounded bands is finite, then q(x) is 
called the finite-gap (g-gap) potential. 

Example 1. We consider the case q(x) = 0, i.e., H = —d 2 /dx 2 . Fix the eigenvalue 
E . Then Hf{x) = Ef(x) is equivalent to f"(x)+Ef(x) = 0. We solve the differential 
equation for dividing into three cases. If E < 0, then we write E = —A 2 and the 
solutions are f(x) = Ae Xx + Be~ Xx for constants A, B, which are unbounded on R 
for {A,B) ytz (0,0). If E = 0, then the solutions are f(x) = A + Bx, which are 
unbounded on M. for B ^ 0. If E > 0, then write E = A 2 and the solutions are 
f{x) = A cos Xx + B sin Ax, which are bouned on IR. Hence we have o~b{H) = (0, 00) 
and o"b(if) = [0, 00). Therefor the potential q(x) = is finite-gap (0-gap). 
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Assume that the potential q(x) is real, smooth, periodic with a period T. Let 
f 1 (x,E), f 2 (x,E) be a basis of solutions to (— d 2 /dx 2 + q(x) — E) f(x) = 0. Then 
fi{x + T, E) and f 2 {x + T, £7) are also solutions and written as 

(2.3) (h(x + T, E) f 2 (x + T, E)) = (fi(x, E) f 2 (x, E))M, 

where M is a 2 x 2 matrix with constant elements. It is known that detM = 1. 
Let t 2 — (trM)t + 1 = be the characteristic polynomial of the monodromy matrix 
M. If |trM| < 2 (resp. |trM| > 2), then E G a b (H) (bounded) (resp. E <£ a b {H) 
(unbounded)). If trM = 2 (resp. trM = —2), then there exists a non-zero periodic 
solution f(x + T) = f(x) (resp. an antiperiodic solution f(x + T) = —f(x)). Hence 
the monodromy caused by the shift of the period (i ^ i+T) implies the boundedness 
or the unboundedness of the solutions to the differential equation. 

The definition of algebro-geometric finite-gap potential is described as follows: 

Definition 2. // there exists an odd-order differential operator 

, , \ 2g+l 2g-l , v 2 g-l-j 

(2.4) A=(±) 



doc J V d/Oc J 



such that 
(2.5) 



a d2 

A -d^ 2+q{x) 



0. 



then q(x) is called the algebro-geometric finite-gap potential. 

Note that the equation [A, —d 2 /dx 2 +q(x)] = is equivalent to the function q(x) be- 
ing a solution to a stationary higher-order KdV equation (see [3]). It was established 
in the 1970s that, under the assumption that q(x) is a periodic, smooth, real function, 
the potential q(x) is finite-gap if and only if q(x) is algebro-geometric finite-gap. 

We now present examples of the (algebro-geometric) finite-gap potentials. Ince [TO] 
established in 1940 that if n e Z> x , u x E R \ {0} and cu 3 G \ {0}, then the 

potential of Lame's operator 

d 2 

(2.6) - — + n(n + l)p(x + uj 3 ), 

dx z 

is finite-gap. In the late 1980s, Treibich and Verdier [22] found that the method 
of finite-gap integration is applicable the elliptic representation of Heun's equation 
for the case lo,li,h,h G Z. Namely, they showed that the potential in Eq. (11.51) is 
an algebro-geometric finite-gap potential if lo,h,h,h G Z. Therefore the potential 
5^i=o ^(^ + + ^i) is called the Treibich- Verdier potential. Treibich and Verdier 

developed the theory of elliptic soliton, Jacobi variety and tangential covering, and 
obtained the results. Subsequently several others [71 [HI [15], [T6], [TJJ US [19] have 
produced more precise statements and concerned results on this subject. 

We obtained further examples of algebro-geometric finite-gap potentials in |20j. 



Theorem 2.1. fl2D]j If M,l ,l h l 2 ,l 3 G Z> 0; <Jj ^ mod 2c^iZ © 2w 3 Z (0 < z < 
3, 1 < j < M), 5j ± 6j> ^ mod 2cjiZ © 2w 3 Z (1 < j < / < M) and Si,..., S M 
satisfy the equations 

3 

(2.7) 2 YjiP'Qi ~ 5 i') + P^i + 6 f)) + + V2) 2 P'(^ + £«;<) = 0, 
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(j — 1, . . . , M), then the potential 

3 M 

(2.8) v{x) = k{k + l)p{x + Ui) + 2 Y,(P( X - + P(x + Sv)), 

i=0 i'=l 

is algebro- geometric finite-gap. 

Note that Eg. (12.71) has appropriate solutions for each M, l ,li,l 2 ,h £ ^>o ( see [20J ) . 
If M = 0, then we recover the result on Heun's equation, and if M = 1, then we 
recover the Treibich's result [2T]. Gesztesy and Weikard [8] developed a theory of 
Picard's potential, and it would be related to our one. 

We introduce a proposition which plays a crutial role of proving Theorem 12.11 
Observe that a product of two solutions to 

(2.9) (j-^ + v{z))f{z) = Ef{z), 
satisfies 

(210) {s? " 4 (t ' w - E) i ~ 3{x - E) = °- 

Proposition 2.2. fl20j,h Under the condition {2. 7] ) ; Eg. li2.1(J\) has a unique non-zero 
doubly-periodic solution E(x,E), which has the expansion 

3 

(2.11) E(x, E) = c (E) + ^Y1 ^ ( E M X + <*) k ~ s 

j=0 jf=0 

M 

+ ^ rf (i ' ) ( J E)(p(x + 5,) + p(x - 5 V )), 
i'=i 

where the coefficients cq{E), b^\E) and d^'\E) are polynomials in E, these polyno- 
mials do not share any common divisors, and the polynomial Cq(E) is monic. 

We set g = deg E c (E). Then the coefficients satisfy deg E bj (E) < g for all i 
and j. Note that the function E(x,E) is frequently used for describing solutions to 
Eq. (12.91) and the monodromy, as we will see in section [3J Write 

9 

(2.12) E(x,E) = Y,a 9 -i(x)E l . 

i=0 

Then ao(x) = 1 and it follows from Eq. (12.10j) that 

(2.13) af{x) - Av(x)d j (x) - 2v\x)a j (x) + 4a' j+x {x) = 0. 
Define the (2g + l)st-order differential operator A by 

d 2 



+ v(x). 



(2.14) A = g {^)|- I (|«,(,))} ^ H=-<L 

It follows from Eq. (12. 131) that [A, H] = 0. Hence the function v(x) is an algebro- 
geometric finite-gap potential, and we obtain Theorem 12.11 Set 

(2.15) Q(E) = E(x, Ef (E - v(x)) + \e{x, E)^^± - \ 
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It is shown by differentiating the right-hand side of Eq. fl2.15l) and applying Eq. fl2.10p 
that Q(E) is independent of x, and Q(E) is a monic polynomial in E of degree 2g + l. 
On the operators H and A, we have the relation A 2 + Q(H) = (see [TTl Proposition 
3.2]). 

Example 2. For the case M = 0, Iq = 2, l\ = l 2 = h — 0, we have 

(2.16) E(x, E) = a (x)E 2 + ai (x)E + a 2 (x) = E 2 + 3p(x)E + 9(p(x) 2 - # 2 /4), 

+ 9{(pW 2 -|)^-pMp'm} 

(2.18) Q(E) = (E 2 - 3g 2 )(E 3 - 9g 2 E/4 - 27<fe/4), 
where g 2 = -4(eie 2 + e 2 e 3 + e 3 ei) and g 3 = 4eie 2 e 3 . 



3. MONODROMY RELATED WITH FINITE-GAP POTENTIAL 

In this section, we review results on solutions to the Schrodinger equation (see 
Eq. (12.91) ) with the algebro-geometric finite-gap potential and its monodromy. Namely, 
we have an integral representation for a solution to Eq. fl2.9p . a monodromy formula 
in terms of hyperelliptic integral, an expression of the Hermite-Krichever Ansatz 
and hyperelliptic-to-ellitic reduction integral formula. The following propositions are 
correct for the potentials in Theorem l2.ll which include the case of Heun's didderential 
equation. 

Proposition 3.1. (Integral representation for eigenfunctions ( \15\ for Heun's equa- 
tion)) Let E(x,E) be the function defined in Proposition \2.2\ and Q(E) be the poly- 
nomial defined in Eq.§2J~B\). Then 



(3.1) 



A(,, £ ) = ^£)exp/^3|M£ 



is a solution to Eq. /\2.9\) . 

Proposition 3.2. (Monodromy formula in terms of hyperelliptic integral ( \17\ for 
Heun's equation)) Let k £ {1,3} ; q k £ {0,1} and E be the eigenvalue such that 
A(x + 2u k , Eq) = (-l)^A(x, E Q ). Then 

(3.2) A(x + 2u k , E) = (-l) q *A{x, E) exp I — - / Jo+£ —dE 



2 



Eq 



-Q(E) 



6 KOUICHI TAKEMURA 

Write 



(3.3) E(x, E) = c(E) + E °f ^ ^ P(* + <* 

i— n n— n \ / 



3 Ji-1 / j N 2j 

EE 

i=0 j=0 

+ J2d {t '\E)(p(x + 5 l ,) + p(x-5 l ,)), 



M 



i'=l 

3 M 

(3.4) a(£)=^a«(i?) + 2^/)(4 

i=0 i'=l 

From Proposition 13.21 we have the following formula, which is expressed as a hyper- 
elliptic integral of second kind: 

(3.5) A(, + a*,^ = (-l)»A(,,£0 •«*(-; f zgMM±j^g) rfjA 

V V-Q(B) / 

for k = 1, 3, where % = £(a>fc) and £(:r) is the Weierstrass zeta function. 

Proposition 3.3. (Hermite-Krichever Ansatz ( \18\ for Heun's equation)) Set 
m , _ s 

(3.6) *(z) = TT(p(*)-pft')). ^(»,«) = ^rx , exp(C(a)g), 

/or z = 0, 1, 2, 3, where o~(x) is the Weierstrass sigma function. There exist polynomi- 
als Pi(E), . . . ,Pq(E) such that, if P 2 (E') ^ 0, then the function A(x,E') is written 
as 

(3 ,) Hx , E ^^M(£p ? ^y Hx , a 

where the values a and k are expressed as 

(3.8) ,M = |ff «a) = |gy=eW, -t|gV=QP0. 

Note that we have the periodicities 

(3.9) A(z + 2u k , E) = exp(-2r/ fc a + 2u k ({a) + 2Kw fc )A(x, -E), 
for k = 1, 3. 

We can obtain hyperelliptic-elliptic reduction formulae by comparing two expres- 
sions (Eq. ( 13.51) and Eq. (13.91) ) of the monodromy. 

Proposition 3.4. (Hyperelliptic-to-ellitic reduction integral formula ( [18] for Heun's 
equation)) Set £ = P 1 (E)/P 2 (E). 
(i) We have 

,3,o, r . =-U B -^E. 

fzij Let cto denote the value of a at E = Eq, where Eq is the value satisfying Q(Eq) = 
0. If oto ^ (mod 2a; 1 Z © 2uj^L), then k = P§{E) ^J—Q(E) / P%{E) is also expressed 
as 
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,3.11) . :, ff> ^ 



Example 3. VFe consider tha case M = 0, Iq = 2, l\ = l 2 = I3 = Recall that the 
functions H(x, £7) and Q(-E) are calculated as 

(3.12) ~(x, E) = E 2 + 3p(x)E + 9(p(x) 2 - a 2 /4), 

(3.13) Q(£) = (E 2 - 3g 2 )(E 3 - 9s 2 £/4 - 27a 3 /4). 

The integral representation for solutions are expressed by substituting them into Eq. $3.1}) . 
The monodromy formula in terms of hyperelliptic integral are written as 

(3.14) A{x + 2u, h ,E) = Hx,E)e*p I 4 f 1 7^ ~ 



2 



392 



-Q(E) 



for k = 1,3. The function A(x,E) admits an expression of the Hermite-Krichever 
Ansatz 

(3.15) A(x, E) = exp (kx) ( bo$o(x, ot) + b\ — <&o(x, a ] 



dx 



and the values a and k satisfy 



(3.16) p(a) = -f- 27 ^, , ' MPZM4-2W4) 
U 9(£ 2 -3 9 ,)' 3Y (B 2 -3 92 ) 

TTie hyperelliptic-to- elliptic reduction integral formula for this case is written as 

1 3£d£ /■« df 

(3.18) 



-(£ 2 - 3a 2 )(^ 3 - f a 2 £ - f a 3 ) ^ V 4 ^ " ^ " ^ 3 
1 /- E (£ 2 -§<? 2 )d£ ft £d£ 



2 



^3 _ 07„_i //'Q/' Z72 



where £ = — (i? — 27g 3 ) / (9(E — 3g 2 )) and k is defined as Eq. A3. 16}) . These formulae 
reduce hyperelliptic integrals of genus two to elliptic integrals. 

4. Finite-gap potential and WKB analysis 

The WKB analysis (the WKB approximation) appears in a semiclassical calculation 
in quantum mechanics, and the WKB analysis is applied for the asymptotic analysis. 
On the WKB analysis of the Schrodinger equation, a solution to the Schrodinger 
equation is obtained as a formal power series in 7]" 1 by introducing a large parameter 
77, and it may have convergent expression by considering the Borel transformation for 
some cases (see [I]). We now consider a solution to the equation 

(4.1) f--f- + V 2 (v(x) - E)) 1>{x) = 0, 



dx 2 

with a large parameter rj. We will find solutions in the form ip(x) = exp(J S(x, rj)dx) 
(4.2) S(x, V ) = S^(x)rj + S (x) + S^x)^ 1 + S 2 (x)rj- 2 + .... 



8 
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We set Q = v(x) — E. Then the function S(x, rf) would satisfy 



(4.3) S(x, V ) 2 + 9 -^- = rjiQ, 



and we have the recursion formula 
(4.4) S- 1 (x) 2 = Q, 2S. 1 (x)S j (x) 



I J2 S ^)Si(x) + 

\k+l=j-l,k,l>0 



dSj-i(x) 
dx 



=3-1 

We set S odd = ^2 j>0 S 2 j-\{x)r] 1 ~^ and S even = £) 7>0 S 2 j(x)r)~ 2: > . Then we have 



(4.5) 



Seven = ~^7^~^~' ^± = ( S °^) eX P f± J S odd dx^ . 



Note that ip± are formal solutions as power series in r]~ . Here we point out that 
Eg. (14.51) resembles to an integral representation of solutions (Eq. ( 13. 1H ) in the theory of 
finite-gap potential. In fact Eq. (13.11) is written as Eq. (I4.5H up to a scalar multiplication 
by setting S odd = y/-Q(E)/E(x, E). Set S = S odd - (dS odd /dx)/(2S odd ). Then an 
equality as Eq. fl4.3|) follows from Eq. (12. 151) . 

We consider the case v(x) = p(x), which is the case of Lame's equation, 

(4.6) (--^ 2 +V 2 (p(x)-E) ) )Hx)=0. 

We have two methods for inventigating Lame's equation, the finite-gap integration 
and the WKB analysis. If rj 2 = 1(1 + 1) for / e Z, then the potential of Eq. (14.61) is 
algebro-geometric finite-gap and we can apply the finite-gap integration for Lame's 
equation. On the limit 77 — > 00, we may apply the WKB analysis. Fedoryuk [5l 
E] investigated Lame's equation, Lame wave equation and Heun's equation by the 
WKB analysis and obtained error estimates for WKB approximation. On the WKB 
expansion of the solutions to Lame's equation, we have 

(4.7) S- 1 (x) = (p(x)-E) 1 / 2 , S (x)- ~"' U) 



A(p(x)-EY 



Sl ( X) =- /: J ! r) L.,- p " {x) 



32{p{x) - Efl 2 S(p(x) — E) 3 / 2 

We consider the leading term of r]" 1 , i.e. S-i(x). Then the monodromy with respect 
to the shift x — > x + 2ujj is written as 

(4.8) S(x + 2ujj) = exp I ± V / (p(x) - E) 1/2 dx + 7]-\. ..) + .. . S{x). 



We will sketch the distribution of the eigenvalues of Lame polynomial for r] — > 00 
and I G Z>o (rj 2 = 1(1 + 1)). Lame polynomial is characterized as a non-zero doubly- 
periodic solution (i.e. f(x + 2uii)/f(x) G {±1} for i — 1,3) to Lame's equation (see 
[231 115]), an d it is essentially a polynomial in z by setting z = p(x). If I is a positive 
integer, then the number of the eigenvalues of the Lame polynomial is 21 + 1 (see 
[23| [2]), and the eigenvalue of the Lame polynomial satisfies rfei < rfE < yfe^, (i.e. 
ei < E < €3) for the case uj\ G M>o and uj^ G \J — HR >0 . If r] 2 = 1(1 + 1), then the 
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eigenvalue E of the Lame polynomial satisfies 

eW 3 + 2tJl 



(4.9) 77 / {p{x) - E) 1/2 dx G ttv^Z, {e 2 <E < d), 

✓>wi+2a;3 

(4.10) 77 / (p(x) - Ef /2 dx G ttV^IZ, (e 3 < E < e 2 ), 



by the asymptotics as rf x — > 0, which follows from the periodicity of the monodromy. 
It is shown under the condition ife\ < rfE < rfet, uj\ G M >0 and cj 3 G \J~^ lM>o that 
a simply-periodic solution is also doubly-periodic. It seems that / pairs of eigenvalues 
merge by the WKB approximation. 

We change a varible in the integral by setting z = p{x). The m-th eigenvalue 
E G (e 2 , ei) of the Lame polynomial from the top should satisfy 



4.11 / x - -dz = n , 

(ei - z){e 2 - z){z - e 3 ) rj 



e3 



and the m-th eigenvalue E G (e 3 , e 2 ) of the Lame polynomial from the bottom should 
satisfy 



, z — E , mn 
4.12 / x - -dz = n . 

(ei - z)(z - e 2 ){z - e 3 ) rj 



as i] 1 — > 0. Hence the number of eigenvalues of the Lame polynomial which is less 
than E is given by 



(4.13) n(E) 



w \J e 3 V ( e l- z )( e 2-z)(z-e :i ) 

as r/^ 1 — > 0. Note that we used the formulae 

1 



-(^-rJi v Z ~ E v s dz), (e 3 <E<e 2 ). 

) , (e 2 < E < ei), 



e2 / E-z 



(414) / ' - 1 = / 

■/e 3 V(e 2 - z)(z - e 3 ) y e: 



y/{e 2 -z)(z- e 3 ) 7 e2 ^/(ei - - e 2 ) 



= it 



^ 3 V(ei - z)(z - e 3 ) / yJea V(ei - z)(z - e 3 ) 
The density of the eigenvalues is written as 

10 ! t - I f 61 / 1 <fe ; (e 3 < -E < e 2 ), 

( 4 - 15 ) = 



f P , 1 ds , (e 2 <E< ei ). 
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Hence we recover the formula given by Borcea and Shapiro [2J by the completely 
different method. We decompose the integral as 

1 f e2 dz 
(4.16) ' 



2 We 3 y/( ei -z)(e 2 -z)(z-e 3 )(E-z) 
1 r 2 dz 



2ny'(e 1 -E)(E-e 3 ) J ea ^(e 2 -z){E-z) 
1 r 2 dz ( 1 I 

+ 7T- 



2vr4 ^fa-z^E-z) \^( ei -z)(z-e 3 ) ^ - E)(E - e 3 ) 
Then the asymptotic of the density of the eigenvalues as E — > e 2 + is written as 
(4.17) 

loo- f 4 ( e 2- e 3A i„„ /" ^ei-ea) ^ \ i / 16(ei-e 2 )(e 2 -e 3 ) 



27r I \/(ei - e 2 )(e 2 - e 3 ) ^/(ei - e 2 )(e 2 - e 3 ) I 27r v /( ei _ e 2 )(e 2 - e 3 ) 
It is also shown that the asymptotic as — > e 2 — is written as 

]„„ /" 16(ei~e 2 )(e 2 -e 3 ; 

(4.18) 1 10g >y (ei-eaXea-i?) 

27r \/( e i - e 2 )(e 2 - e 3 ) 
Hence the density of the eigenvalues as E — > e 2 has logamithtic singularity. If e 2 = e 3 , 
then the density is written as l/(2^(e 2 — e±)(E — e 2 )) and the type of the singularity 
at E = e 2 changes. 

We have another way to insert a large parameter, we regard the eigenvalue E as 
a large parameter while fixing the potential. Namely, we consider the Schrodinger 
equation 

(4.19) (-^ + (v{x) - E)^j iP{x) = 0, 
by setting E = —rf, 

(4.20) ip(x) = exp(r)x + ip^x)^" 1 + il) 2 (x)ri~ 2 + 4'3( x ) r l~ 3 + V^C^) 7 ? -4 + •••)• 

and we regard 7] as a large parameter while fixing the potential. We now restrict to 
the case of Lame's equation v(x) = 1(1+1) p(x). Then the coefficients ipi(x), ilfy(x), .. . 
are expressed as 

(4.21) Mx ) = Jl±R C{x)i Mx) = - l l±H p ( x ), 

l 2 (l + l) 2 ( l 2 (l + l) 2 l(l + l)\ „ , 

ibJx) = i >—g X + 1 '— + _V L p'( X ) 

, . . l 2 (l + l) 2 (l 2 {l + l) 2 l(l + l)\ „. \ 

ibJx) = — —g 2 + — p (x), .... 

v V ; 96 y V 48 16 J 

Hence we have the asymptotics of the monodromy 

(4.22) if}(x + 2ui) = ip(x) exp (2^77 - 1(1 + l)ViV' X ~ — ^f^ V~ 3 + ■■■ 

V 48 

where rji = C( w i) (i — 1 ; 3). Grosset and Veselov [9] studied the asymptotics of the 
monodromy by applying the elliptic Faulhaber polynomials, which is related to the 
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soliton theory of the KdV equation. They considered the asymptotics of the densities 
of states for Lame's equation, and applied it for computation of the polynomials 
—2rjka(E) + 2u)kc(E) in Eg. ( 13.5(1 . Eq.f l4.22p is also obtained as a consequence of the 
paper [9]. 

We expect further results among the finite-gap potential, WKB analysis and Heun's 
(Lame's) equation. 
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